Polariton dynamics and Bose-Einstein condensation in semiconductor microcavities 
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We present a theoretical model that allows us to describe the polariton dynamics in a semi- 
conductor microcavity at large densities, for the case of non-resonant excitation. Exciton-polariton 
scattering from a thermalized exciton reservoir is identified as the main mechanism for relaxation 
into the lower polariton states. A maximum in the polariton distribution that shifts towards lower 
energies with increasing pump-power or temperature is shown, in agreement with recent experi- 
ments. Above a critical pump-power, macroscopic occupancies (5 x 10*) can be achieved in the 
lowest energy polariton state. Our model predicts the possibility of Bose-Einstein Condensation of 
polaritons, driven by exciton-polariton interaction, at densities well below the saturation density for 
CdTe microcavities. 

PACS numbers: 71.35.-)-z 
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I. INTRODUCTION 



Polaritons are quasiparticles created by the strong 
coupling between excitons and photons, and behave as 
composite bosons at small enough densities. In the 
last decade a huge experimental and theoretical effort 
has been dedicated to the search for quantum degener- 
acy effects in microcavity polaritons. They were soon 
proposed as candidates for the formation of a Bose- 
Einstein condensate, due to their small density of states 
(DOS) However, a bottleneck effect due to the slow 
polariton-acoustical phonon scattering was predicted |^] 
and observed in angle-resolved experiments Stim- 
ulated scattering due to the bosonic nature of polari- 
tons has been demonstrated Q in pump-probe experi- 
ments. Parametric amplification and parametric oscilla- 
tion, under resonant excitation, has been observed [^||^] 
and was well explained by coherent polariton-polariton 
scattering In GaAs microcavities, experiments 

have shown the role of the polariton-polariton scattering 
in the nonlinear emission after non-resonant excitation, 
either in continous wave |lo|-[l^ , or time-resolved exper- 
iments . The growth of II- VI microcavities leads to 
new exciting possibilities due to the larger stability of 
the exciton. The study of these samples has allowed 
for the observation of stimulated scattering in the strong 
coupling regime under non-resonant excitation conditions 
iQ . The robustness of the polariton in these samples has 
been recently shown by the demonstration of parametric 
amplification at high temperatures p5| . 

The theoretical study of semiconductor microcavities 
has followed two different lines. A fermionic formalism, 
including carrier-carrier correlation, explained the lasing 
by population inversion at densities above the saturation 
density, where the coupling between the exciton and the 
photon disappears (weak couphng regime) |]l^,^. On 



the other hand, a bosonic picture, including polariton- 
polariton interaction, has been developed. It must be 
stressed that recent experiments agree with the predic- 
tions made by the bosonic formalism (bottleneck effect, 
stimulated scattering) under excitation conditions such 
that the density of polaritons remains below the satura- 
tion density. The polariton-polariton interaction, given 
by the usual Coulomb exchange interaction between ex- 
citons, and the saturation term, has also explained quan- 
titatively many recent observations |^-|^. 

In this paper we present a theoretical study of the evo- 
lution of the polariton population at large densities with 
a bosonic description. We use the semiclassical Boltz- 
mann equation, developed by Tassone et al. Q, as the 
starting point for the description of the polariton dy- 
namics. Exciton-polariton (X-P) scattering is identified 
as the most important mechanism for relaxation towards 
the lower polariton branch. We deduce a simplified model 
in which the high energy excitons are considered as a 
thermalized reservoir. This simplification allows us to 
describe all the microscopic lower polariton states, and 
the evolution of the polariton population towards Bose- 
Einstein Condensation (BEC) |T9| ]. 

X-P scattering is shown to produce a polariton pop- 
ulation with an occupancy that has a maximum at a 
given energy below the exciton energy. The energy at 
which the polariton distribution peaks is shown to de- 
crease linearly as pump-power or sample temperature 
is increased, in such a way that the dip in the polari- 
ton dispersion behaves as a polariton trap in phase-space 
pO{ . This conclusion is valid for all semiconductor mi- 
crocavities. In particular, the same behavior has been 
observed in a recent experiment in a GaAs microcavity 
pO| . However, in GaAs, densities comparable to the sat- 
uration density are created, as one increases the pump- 
power, before the threshold for BEC is reached. Thus, for 
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a clean observation of all the phenomenology described 
in this paper, better material parameters than those of 
GaAs microcavities are desirable. We predict that in 
CdTe microcavities, where excitons are more stable due 
to their smaller Bohr radius, the threshold for BEG can 
be reached for densities well below the saturation den- 
sity. Thus, we choose, for concreteness, CdTe parameters 
for the detailed calculations shown below. Our model 
shows that the relaxation towards the ground state is 
only possible for densities comparable to the saturation 
density in GaAs microcavities. In order to achieve BEG 
with smaller polariton densities, an additional scattering 
channel that those considered here should be included. 
Electron-polariton scattering in doped microcavities has 
been recently proposed (2l| as an efficient mechanism 
that allows BEG for such small polariton densities. 

The paper is organized as follows. In section 11 we 
present the theoretical framework. A set of differential 
equations is obtained for the description of the micro- 
scopic lower polariton levels, and the exciton density and 
exciton temperature. In section III, we apply this model 
to the study of the dynamics of the polariton popula- 
tion for different pump-powers and temperatures, after 
non-resonant, continuous excitation. The case of Bose- 
Einstein Gondensation into the lowest energy polariton 
mode is considered in section IV. Finally, in section V, 
we list our main conclusions, and discuss the comparison 
with experiments. 



II. RATE EQUATION SIMPLIFIED MODEL 

A. Simplified model for the exciton-polariton 
scattering 

At densities such that na^ << 1, where qb is the 
2-D exciton Bohr radius, excitons can be described as 
interacting bosons. From typical parameters for GdTe 
QW's (exciton binding energy Eb « 25 meV, and elec- 
tron and hole effective masses me = 0.096, nih — 0.2), 
we get a 2-D exciton Bohr radius as — ^7 A, and 
— 4.5 X 10^^cm~^. On the other hand, at large ex- 
citon densities, carrier exchange interactions and phase- 
space filling effects destroy the coupling between excitons 
and photons. The density at which the exciton oscillator 
strength is completely screened is the saturation density, 
Usat- Following Schmitt-Rink et al. [ |2^ , we can estimate, 
for a GdTe microcavity. 



risat 



= 6.7 ■ lO^'cm-^ 



ira 



(1) 



where a^^^ — as/^. In InGaAs and GaAs QW's, Ugat 
would result to be 6.6 x 10^°cto^^ and 1.3 x lO^^cm"^, 
respectively. Along this paper we consider the case of 
GdTe in a range of polariton densities well below Usat, 



so that the microcavity remains in the strong coupling 
regime. 

In the strong coupling regime, the Hamiltonian in- 
cludes the coupling between excitons and photons 
through the polariton splitting, Qp p3[: 



(2) 



^k' ^k ('^k' ^k) creation operator and dispersion re- 

lation for the bare exciton (photon) mode. The residual 
Goulomb interaction between excitons, and the satura- 
tion of the exciton oscillator strength, can be described 
by the effective exciton-exciton, and exciton-photon in- 
teraction Hamiltonian ]2^ : 

M. 



^1= J2 ^^L^L^ki+q&k^-q 



ki.k2 



-^^ki^L^ki+qOka-q + Zi-C. (3) 

where S represents the quantization area. We consider 
here the small momentum limit for Hj, an approximation 
that is valid for the range of small wave-vectors consid- 
ered below. The matrix elements in Hj are given by: 



M^^ « 6 Ssal 
Csat ~ 1-8 il.pa% 



Hq can be diagonalized in the polariton basis: 



Pip,k = &k + C^ak 



Pup,k — AT^^'&k -r 



(4) 



(5) 



where X^^^^, C^^^^ are the Hopfield coefficients, that 
represent the excitonic and photonic weights in the po- 
lariton wave- function for the lower/upper branch, respec- 
tively ||2^ . By means of the transformation (^) we can ex- 
press Hj in terms of polariton operators. We will neglect 
the upper polariton branch, an approximation that will 
be justified later. We express the operators 5k, Ck in (||) 
in terms of p/p,k, Pup,k, and retain the terms correspond- 
ing to interactions within the lower polariton branch only. 
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25'^ki,kr,k3,k4^'/p,ki^'/p,k2?''p,k3^'ip,k4' i^) 



where 
1, 



ki,k2 
k3,k4 



lp~lp 



2 ^1 ,k2 ,k3,k4 

^ M X^P V-P y'-p 
2 



M.. x;^x;^x-:^x^^ 



f^ip Y^P y'p v'p 



^ v^p y'-p r"-p y 



5ki+k2,k3-|-k4 



(7) 
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is the dispersion relation for the lower polariton 
branch. It is plotted, assuming typical parameters for 
a semiconductor microcavity, in Fig. |l[ For energies 
higher than the bare exciton energy, it matches the bare 
exciton dispersion. 

The theoretical model developed by Tassone et al. 
p8[ describes the evolution of the polariton popula- 
tion by means of a semiclassical Boltzmann equation 
that includes polariton-acoustical phonon and polariton- 
polariton scattering. A grid uniform in energy is consid- 
ered for the description of the polaritonic levels. Here, 
we are interested in the evolution of the system towards 
a Bose-Einstein condensate. It is well known that in 2- 
D systems, EEC is possible at finite temperatures, only 
for systems with a finite quantization length, Lc- As a 
result, the critical density and temperature depend log- 
arithmically on this length scale. Due to the peculiarity 
introduced by the 2-D character of semiconductor mi- 
crocavities, the grid uniform in energy considered in fl^ ] 
is not well suited for the study of the transition towards 
EEC. Instead, a uniform grid in k-space, related with the 
inverse of the quantization length, , has to be consid- 
ered. Unfortunately, the number of levels that one has 
to consider for the description of the polariton dynam- 
ics is dramatically increased, and a numerical calculation 
that includes all the possible polariton-polariton scatter- 
ing processes is not possible. Thus, some simplification 
has to be done, which still allows to describe satisfacto- 
rily the polariton dynamics, and to predict the possibility 
of Eose-Einstein Condensation. 

In order to develop a simplified model for the case of 
non-resonant excitation, we make the following assump- 
tions: 

(1) The exciton population above the bare exciton en- 
ergy can be considered as a reservoir at a given temper- 
ature, that follows a Maxwell-Eoltzmann distribution. 

(2) The main scattering mechanism at large densi- 
ties for relaxation into the lower polariton modes is the 
one shown in Fig. ^, in which two excitons scatter in 
the reservoir, and the final states are a lower polariton, 
and another reservoir exciton. We will call this process 
exciton-polariton scattering. 

A scattering process in which two excitons are initially 
in the exciton reservoir, and a lower and upper polariton 
are the final states is also possible, but has a much smaller 
probability. This is due to the fact that in such a colli- 
sion, an upper polariton, instead of an exciton, plays the 
role of a final state. Thus, the scattering rate is reduced 
by the ratio between the upper polariton and bare exci- 
ton DOS, Pupl Px ~ 4 X 10^^. In fact, for negative or zero 
detunings, the upper polariton population remains negli- 
gible, as shown experimentally jl4| and theoretically p8| . 
The scattering of excitons with acoustical phonons going 
into lower polariton states is also much slower than the 
X-P scattering already at exciton densities of lO^cm"^, 
and the lattice temperatures considered here. 



The phase-space is, thus, divided into two pieces. 

• The exciton reservoir, which describes the excitons 
directly created by the non-resonant pump. We 
take as zero energy the bare exciton energy, so that 
the exciton reservoir includes all the levels with 
e > 0. These levels are weakly coupled to the pho- 
tonic modes, so that they can be approximated as 
bare excitons. The occupation numbers are given 

by, 



Nf = iV^e" 



is given by 



2'Knx 
p-xknTx ' 



(8) 



(9) 



where rix is the exciton density, px is the bare ex- 
citon density of states, 



Px = {me + mh)/h^ 



(10) 



and Tx is the temperature of the exciton reservoir. 
The assumption of a rigid thermal exciton distri- 
bution is well justified, as discussed below, for lev- 
els above e = 0, due to the fast exciton-cxciton 
scattering that thermalize the exciton population 
at large densities [ p^j26| . On the other hand, 
a Maxwell-Eoltzmann distribution, rather than a 
Eose-Einstein distribution can be used to describe 
the exciton reservoir, provided that N'^ << 1. 

« 0.2 for the largest densities considered here, 
so that quantum degeneracy effects are negligible. 
The main mechanism for radiative losses in the ex- 
citon reservoir is due to the coupling to the leaky 
modes, that we include by considering an exciton 
hfetime, = lOOps 

• The lower polariton branch (LPB), where the en- 
ergy levels are below the bare exciton energy (e|f < 
0). Exciton and photon levels in the LPE are 
strongly coupled, and radiative losses are much 
more important, due to the photonic weight of 
these states. Thus, we expect a nonequilibrium po- 
lariton distribution that depends strongly on tem- 
perature and density. On the other hand, the 
abrupt dispersion relation implies a very small ef- 
fective mass for the LPE states, mip « 4 x lO^^rn^, 
leading to the possibility of Bose-Einstein Conden- 
sation. 

With this separation of the phase-space we are neglect- 
ing the details of the polariton population near e « 0, 
that is, at the knee in the polariton dispersion relation, 
where the LPE and the quadratic dispersion of the bare 
excitons merge together (see Fig. |l|). The contribution 
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from these states is negligible already at the lowest tem- 
peratures considered here {4:K, ksT ~ 0.35meT/), be- 
cause the main part of the exciton distribution is in the 
weak coupling region. 




FIG. 1. Exciton-polariton scattering process from a ther- 
malized exciton reservoir into the lower polariton branch. 

We describe the LPB states by a grid in k-space. The 
exact quantization area is not well determined in exper- 
iments, due to inhomogcncities in the spot, and surface 
defects. The relevant quantity is the quantization length, 
that gives the order of magnitude for the distance be- 
tween the energy levels. We consider first Lc — bOfim 
(a typical excitation spot diameter), and plane waves for 
the microscopic polariton states, so that the quantiza- 
tion area is simply S — L^. The allowed wave- vectors 
are given by 

277 

k'P^{n^^ + nyy) — , (11) 

where n^^Uy are integers. A similar method has been 
employed to describe the microscopic polariton levels in 
a microcavity etched into a microscopic post structure 
(Lc = 2/iTO) In that case, the LPB is reduced to 

a single microscopic level. On the contrary, in our case, 
Lc is larger, and the number of lower polariton levels re- 
sults to be of the order of 10''. The main results of our 
work do not depend strongly on the exact value of Lc, 
but only on its order of magnitude, as we show in the 
last section. The weak dependence of our results on Lc 
can be easily understood from the peculiarities of EEC 
in 2-D (see Appendix B). 

We write now a set of rate equations for iV^^ (occupa- 
tion numbers in the lower polariton branch) and n,^. (the 
density of the thermalized exciton reservoir). We con- 
sider the injection of polaritons from isotropically dis- 
tributed excitons, so that the occupation numbers N^^ 
will depend on the absolute value of the in-plane wave- 
vector only. The rates for scattering from the thermalized 



reservoir into and out of the LPB, due to the interaction 
Hamiltonian (||) are calculated in Appendix A. They lead 
to a simplified version of the semiclassical Boltzmann 
equation: 

dN'P _ 
dt ~ 

drix 
~dt ~ 

-^xTlx+Px- (12) 

P|f , Tx are the radiative losses in the LPB and exciton 
reservoir, respectively: 

Tc 

= -, (13) 

Tx 

where Tc is the lifetime for the bare photon within the 
cavity. We consider = Ips, a typical value for high 
quality microcavities. Px represents the non-resonant 
pump that injects excitons directly into the exciton reser- 
voir. If relaxation into the LPB is slow enough, the life- 
time for an exciton in the exciton reservoir is given by 
Tx, so that a given pump px = Ux/tx = Ux/lOOps cre- 
ates an exciton density Ux- If scattering into the LPB is 
fast enough, however, the density created by the pump is 
smaller, as we will see below, dg^^ is the degeneracy for 
the polariton level corresponding to the wave-vector fc, 
and is calculated according to the distribution obtained 
from Eq. (|l]). 

W^'^^Ux are the rates for X-P scattering into 
and out of the LPB. They are given by the expressions 
(see Appendix A): 

= TMipxe^<"^''^^^, (14) 
n 

where 

Mk ^ MxxX'P + a.atCf . (15) 

The dependence of the rates on temperature can be 
qualitatively understood by the phase-space restrictions 
for relaxation into and out of the LPB (Fig. Let 
us consider two excitons in the exciton reservoir with 
a wave- vector kf and energy ef, that scatter to a LPB 
state with energy and to a high energy state with k^ , 
e^. Taking into account that the wave- vector of the low- 
energy polariton is negligible, when compared with the 
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wave-vectors in the reservoir, simultaneous energy and 
momentum conservation implies that: 



2 ef 



2 \K\ = Wf\ 



4 6?= = 



ipi 



(16) 



Eq. ( p^ ) means that two excitons need an energy lejf |/2 
to be the initial states of a process in which one of them 
falls into the LPB. On the other hand, an exciton needs 
an energy 2 \e^\ in order to scatter with a lower polariton 
and take it out of the LPB. 

The polariton distribution in the steady-state regime 
will be given by the set of equations (|2|) with dNjF /dt — 
0, drix/dt = 0. In the steady-state, the occupation num- 
bers in the LPB are given by: 



-dp 



(17) 



If we neglect the dependences on k of the matrix element 
Mfc and the radiative losses F^?', then the maximum in 
the polariton distribution has a peak at the energy: 
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Eq. (g^ implies an approximately linear relation be- 
tween the exciton reservoir temperature and the energy 
position at which the system prefers to relax into the 
LPB (that is, the maximum energy loss in the relaxation 
process). 



B. Model for the evolution of the exciton reservoir 
temperature 

The strong dependence of the X-P scattering into the 
LPB on Tx implies that, in order to complete our model, 
we have to describe properly the evolution of the exci- 
ton reservoir temperature. During the relaxation process 
into the LPB, shown in Fig. |l|, a high energy exciton is 
injected into the exciton reservoir. This way, during the 
relaxation process, the exciton reservoir is heated. On 
the other hand, the scattering with acoustical phonons, 
tries to keep the exciton reservoir at the lattice temper- 
ature. In order to describe the evolution of the exciton 
reservoir temperature, a new variable has to be added 
to this set of equations. We introduce e^, the density 
of energy for the exciton reservoir, that is given, for a 
Maxwell-Boltzmann distribution, by the expression: 



s 



(19) 



The evolution of the energy allows us to obtain the ex- 
citon temperature, at each step in the evolution of the 
polariton population. 



However, we still need to make an assumption about 
the exciton distribution that the pump directly injects 
into the reservoir. Non-resonant experiments usually cre- 
ate a population of high energy electron-hole pairs, that 
relax towards low energy states, after the process of ex- 
citon formation 29|. It must be stressed that for the 
case of CdTe QW's this process is particularly fast, due 
to the possibility of exciton formation through the emis- 
sion of LO phonons, whose energy is 21 meV, below the 
exciton binding energy (25 meV). The time scale for ther- 
malization due to exciton-exciton scattering is very fast 
(« Ips), as compared with the time scale for relaxation 
into the LPB (w 20ps, for the larger exciton densities 
considered here). Thus, it is very well justified to de- 
scribe the non-resonant pump, as a source that injects 
thcrmalizcd excitons at a given pump-temperature into 
the exciton reservoir. We still need to specify which 
is the initial pump-temperature. Due to the large un- 
certainties in non-resonant experiments it is very dif- 
ficult to make predictive calculations about this initial 
condition. We choose instead to make the simplest as- 
sumption: we consider that the pump injects excitons at 
the lattice temperature, T^. For large exciton densities 
(18) (above 5 x lO^^cm"^) the processes of reservoir heating 
by X-P scattering and cooling through exciton-phonon 
scattering, govern the final exciton temperature, and our 
results do not depend on our assumption for the pump- 
temperature. For small temperatures (below 10 K), and 
densities, our results depend on this assumption. In a 
real experiment, the exciton temperature can be larger 
than the one considered here, if the cooling of the initial 
non-resonantly created excitons is not fast enough. 

It must be pointed out that a very efficient experi- 
mental setup for non-resonant excitation of a cold exci- 
ton distribution in a semiconductor microcavity has been 
demonstrated by Savvidis et al. [|o|. In this experiment, 
non-resonant excitation has been achieved with small ex- 
cess energies (below the stop band), due to the finite 
transmittance of the Bragg refiectors. The experimental 
results presented in are in agreement with a model in 
which excitons relax from a thcrmalizcd distribution into 
the lowest energy states, and justify all the assumptions 
made above. 

We have to quantify the three contributions for the 
evolution of the exciton temperature listed above: 



(i) Heating by the scattering to high energy exciton lev- 
els. 

The total energy is conserved in the X-P scattering to- 
wards the LPB. Thus, the energy lost by the polaritons 
falling into the lower energy states must be gained by the 
exciton reservoir. The following expression gives the rate 
of change of Cx by each X-P scattering process: 



dex 
~dt 



k 
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X {Wrnlil + iV^) - Wrn.N[^). (20) 

The heating of the exciton gas in the relaxation to- 
wards the LPB was also shown theoretically in for 
the case of a microcavity post. 

(a) Cooling produced by acoustical phonons. 

The exciton-phonon scattering tries to keep excitons 
at the lattice temperature, T^, against the heating con- 
sidered above We divide the exciton Maxwell- 
Boltzmann distribution into levels equally spaced in en- 
ergy with a given AE << kBT^^kBTx- The probability 
for absorption or emission of acoustical phonons gives 
a scattering rate W^!^^{Tl) for the transition between 
states with energies ei,ej- These rates depend on the 
phonon distribution through T^, and can be evaluated 
from the exciton-acoustical phonon interaction Hamilto- 
nian. The evolution of the occupation numbers, N^, in 
the exciton reservoir leads to an expression for the rate 
of change of the total energy: 



de, , _ Px^Taf .^'^^^ , 



dt 



\ph-- 



I 3 

AE n. 



ksTx 



E 
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(21) 

For the calculation of W''^ we have used the exciton- 
acoustical phonon coupling through the deformation po- 
tential interaction. These rates have been presented in 
several previous works and we will not re- 

write them here. We consider deformation potentials 
fle = 2.161^, at = 5. lei/ for CdTe [|o), and an 80A QW. 
In Eq. (pl|), stimulated terms are neglected, due to the 
small occupation numbers in the exciton reservoir. The 
rate of change for Cx results to be linear in n^^ and de- 
pendent on both Tl and T^. 

(Hi) Cooling to the pump temperature. 

The injection of excitons into the reservoir by a pump 
that is considered to be at the lattice temperature tries 
to keep the exciton reservoir at T^. The terms that cor- 
respond to pumping and radiative losses in Eqs. (p^), 

\pu — Px ^x^x^ (22) 



dt 



imply that the rate for the evolution of the occupation 
numbers has the contribution: 



dt 



2Trpx 



\pu- 



RxknTh 



^ 27rn:r _ 
1 X — ; e 

PxKbTx 



where the condition that the pump injects excitons at 
into a Maxwell-Boltzmann distribution at Tx is included. 



The rates in Eq. (23) leads to a new contribution to the 
evolution of Cx- 



dt 



k 



,dNl 
dt 



PxksTL 



-kBTx. (24) 



Note that this contribution cools the exciton reservoir if 
Tx > Tl. 



The set of first-order differential equations ([l^), to- 
gether with ( ^ , 24), leads to a closed description 



for the dynamics of the exciton reservoir and lower po- 
laritons. We have integrated numerically these equa- 
tions until the stationary solution {dN''P /dt ~ dux/dt — 
dcx/dt — 0) is reached. The obtained solution describes 
the steady-state for the case of a continuous non-resonant 
pump. 



(23) 



C. Comparison between our simplified model and 
the complete semiclassical Boltzmann equation. 



We have checked whether the polariton distribution 
obtained with this simplified model agree with the one 
obtained in a detailed calculation performed with a 
Boltzmann equation in which the simplifying assump- 
tion of a Maxwell-Boltzmann distribution for the lev- 
els in the exciton reservoir is not made. In Fig. |^ we 
present the comparison between a complete calculation 
that follows reference |l^ , with a grid uniform in energy 
{AE = O.lmeV), in which all the polariton-polariton and 
polariton-phonon scattering processes are included, and 
our simplified model. We expect that the results from 
both calculations coincide under conditions for which 
Bose-Einstein Condensation (that is, the macroscopic oc- 
cupancy of a single microscopic mode) is not possible. 
On the contrary, in the range of densities near or above 
BEC, comparison between both models would not be al- 
lowed because the system would not be free to relax into 
a single microscopic mode. Thus, a continuum density of 
states would not be able to describe properly the LPB, 
and a grid uniform in energy could not be applied. In 
Fig. ^ a few pump-powers have been chosen, for which 
the comparison is possible. It can be clearly seen that 
both calculations yield identical results, so that, we can 
conclude that our simplified model captures the essen- 
tial physics in the relaxation from the exciton reservoir 
towards the lower energy polariton states. 
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-2 
Energy (meV) 

FIG. 2. Polariton distribution calculated by means of the 
complete semiclassical Boltzmann equation, as in Ref. [18] 
(points), and polariton distribution calculated by means of 
our simplified model (continuous line). We have considered 
a CdTe microcavity with zero detuning and Tl ~ lOK. The 
curves, from bottom to top, refer to pump-powers from 1 to 
5, in units of lO^^cm-VlOOps. 



III. EVOLUTION OF THE POLARITON 
DISTRIBUTION WITH PUMP-POWER AND 
TEMPERATURE 



with Tx, smaller occupation numbers are predicted for 
higher temperatures. At high T^, the bottleneck is sup- 
pressed and the polariton distribution becomes almost 
flat. The difference between high and low was already 



observed in 1 11 1 , where a flat distribution was evidenced 
for pump-powers where the X-P scattering is the domi- 
nant relaxation mechanism. 
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FIG. 3. Evolution of the polariton distribution for 
Tl ~ 4, 10, 15 . . . SOTf, flp = 10 meV, zero detuning, and 
VlOOps. 



A. Evolution with temperature 

In the previous section, we obtained an estimation for 
the evolution of the polariton population as the exciton 
temperature is increased. A rough linear relation be- 
tween the position of the maximum in the distribution 
and Tx was predicted. Now we present a complete calcu- 
lation with our model that confirms this estimate. This 
effect can be observed for a different range of polariton 
splittings and detunings. However, the main parameter 
that governs the relaxation from the exciton reservoir is 
the energy difference between the bare exciton and the 
bottom of the LPB. This energy distance can be inter- 
preted as the depth of the polariton trap formed by the 
dip in the polariton dispersion. Within this paper we will 
consider, for concreteness, the case of zero detuning (that 
is, E,j: (bare exciton energy) — (bare photon energy)), 
so that the lowest polariton mode has an energy Q.p/2 
below Ex- 

In Fig. ^ we show the evolution of the polari- 
ton distribution for Vlp = \QmeV . We choose — 
lO^'^cm^^/lOOps. Under such excitation conditions, X-P 
scattering is slower than radiative recombination in the 
exciton reservoir, and Ux ~ TxPx- The heating of the ex- 
citon gas is also negligible, so that ~Tl. The peak in 
the polariton distribution moves towards lower energies 
as temperature is increased. Due to the fact that the rate 
given in Eq. (jl^) for scattering into the LPB decreases 



The evolution of the peak in the polariton distribution 
is shown in Fig. ^ for a range of polariton splittings. The 
linear dependence predicted by the estimation in Eq. ( p^ ) 
is clearly observed. For high enough temperatures and 
small energy distance between the lowest polariton en- 
ergy and the exciton reservoir, the maximum in the po- 
lariton distribution can reach the bottom of the LPB. 
However, the occupation numbers are still smaller than 
one, and the situation is not able to produce BEC. 
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FIG. 4. Dependence of the energy of the peak in the po- 
lariton distribution on Tl for px ~ lO^^cm'^/lOOps, and dif- 
ferent polariton splittings. 
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Due to the isotropy of the lower polariton distribution, 
the emission at a given energy ej^^^^. , corresponds to emis- 
sion at a given angle with the growth axis. The movement 
of e^^ax towards lower energies implies that the sample 
emits light preferentially in a ring, at a given angle that 
decreases with temperature. The recent experiment re- 
ported by Savvidis et al. |l^] has observed the ring emis- 
sion with the linear dependence on lattice temperature 
described above, and shows that the X-P scattering from 
a thermalized exciton reservoir considered here, is the 
main scattering mechanism for relaxation into the LPB. 



B. Evolution with Pump-Power 

As shown in Eq. (p^, ej^^^. depends logarithmically 
on the exciton density, and linearly on Tx- At large den- 
sities, however, the heating of the exciton reservoir must 
be properly taken into account. When occupation num- 
bers N^^ > 1 are achieved, the X-P scattering mechanism 
shown in Fig. |l| is stimulated, and the injection of high 
energy excitons becomes more and more important. In 
this situation, the heating of the exciton reservoir gov- 
erns the evolution of ej^^^, as Px is increased, because 
for higher , the polariton distribution has a maximum 
at lower energies. This way, the heating of the exciton 
reservoir allows the excitons to fall into the polariton trap 
formed by the dip in the polariton dispersion. 

Fig. H shows the evolution of the polariton distribu- 
tion as one increases the pump-power. The distribution 
can reach the lowest energy state and large occupation 
numbers can now be achieved. For large enough powers, 
the lowest energy state can have the main part of the po- 
lariton population, a situation that will be described in 
detail in section IV. For intermediate pump-powers, oc- 
cupation numbers larger than one are obtained for modes 
with k ^ 0. Thus, our model predicts the stimulation to- 
wards these modes before BEC is achieved. 
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FIG. 5. Evolution of the polariton distribution for 
Tl = WK, flp = lOmeV, and = 1,2,5,8, 

15 X 10^°cm"^/100ps (from bottom to top). 



In Fig. ^ we show the evolution of e^^axi one 
increases Px, for different flp's. Three stages can be 
clearly distinguished. From px ~ lO^cm^^/lOOps to 
Px » 10"'^''cm~^/100ps, the exciton density created by the 
pump is small enough for the heating of the exciton reser- 
voir to be neglected. Under such conditions, the evolu- 
tion of e^j^ax is governed by the logarithmic dependence on 
Tlx ~ PxTx (see Eq. |l^)). Above Px ~ 10^°CTO~^/100ps, 
the heating of the exciton gas is important, and the evo- 
lution of ej^jjj. is governed by the evolution of Tx , through 
the linear dependence shown in Eq. (|l^). In this regime, 
^max results to depend linearly on px, until it saturates 
to the bottom of the LPB. 



15 




FIG. 6. Dependence of the energy of the peak in the polari- 
ton distribution on pump-power, for Tl ~ WK, and different 
polariton splittings. 



The heating of the exciton reservoir is shown in Fig. 
. For large pump-powers excitons are no longer at the 
temperature T^, at which they are initially injected. The 
raise of Tx becomes more important for larger flp due 
to the highest energies of the excitons injected into the 
reservoir by the X-P scattering (see Fig. |l|). 
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FIG. 7. Dependence of the temperature of the exciton 
reservoir on pump-power for Tl = lOK, and different po- 
lariton splittings. 
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FIG. 8. Dependence of the total exciton reservoir and lower 
polariton density, n^p, on pump-power, for Tl ~ lOK. The 
results for flp = 6, and lOmeV are indistinguishable. 



At large px, X-P scattering is faster than the radia- 
tive recombination given by t^, an it becomes the main 
mechanism for losses in the exciton reservoir. Thus, the 
exciton-polariton density, Uxp (that is, the density of ex- 
citons in the reservoir plus the density of lower polari- 
tons), created by the pump is smaller than PxT-x- Fig. ^ 
shows the evolution of the steady-state polariton density 
as Px is increased. For all the three polariton splittings 
considered the evolution of the exciton-polariton density 
is very similar. 



The shift of e^ax towards lower energies results also 
in the emission in a ring at an angle that decreases as 
pump-power is increased. In the experiments performed 
by Savvidis et al. this effect was observed, and a lin- 
ear relation between e\^^ax pump-power was reported. 
In this experiment, the extrapolation of the linear de- 
pendence of ej^ttj. on Px towards small powers shows an 
energy offset, in agreement with the theoretical result 
(see Fig. ^). Angle-resolved measurements allow, thus, 
to monitor the evolution of the polariton population to- 
wards the ground state. 

As shown in Fig. |^ for intermediate pump-powers the 
system does not relax down to the ground state, but 
stimulation (occupancy larger than one) is achieved in 
the scattering into intermediate modes with k ^ {]. The 
stimulated modes will thus emit light at angles 6 ^ Q. 
Fig. ^ presents the evolution of the occupation numbers 
for the modes corresponding to different angles, as pump- 
power is increased. The evolution of the occupancy at 
each angles matches the experimental results in pO[. 



Uxp at which 



Fig. H and Fig. || allow to estimate the 
the relaxation into the ground state is possible. In Fig. 
H pump-powers above 5 x 10^°cm~^/100ps correspond 
to polariton densities larger than 4 x 10^°cm^^. Thus, 
in GaAs microcavities, the weak coupling regime will be 
reached before ej^^^. reaches the bottom of the LPB dis- 
persion, due to the small n^at, and BEG is not possible. 
For the clean observation of the complete evolution of the 
polariton distribution towards the ground state, CdTe 
microcavities are promising candidates because Uxp re- 
mains well below Usati before the system relaxes towards 
the ground state. 
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FIG. 9. Evolution of the occupancy of the modes corre- 
sponding to emission angles = 0, 5 ... 25 degrees (from 
right to left), as pump-power is increased. Tl ~ lOK, 
= WmeV. 

A quadratic relation is predicted, followed by an ex- 
ponential growth that saturates at large powers. Stimu- 
lation can be observed, at polariton densities below the 
BEC, at given angles, which correspond to the modes 
where the polariton distribution peaks, as shown in Fig. 
^. The saturation of the emission at a given angle is re- 
lated to the evolution of the polariton population that 
" shifts downwards" in energy, due to the heating of the 
exciton reservoir. Thus, the emission at a given angle sat- 
urates when the peak in the polariton distribution turns 
to smaller angles, as pump-power is increased. 



IV. BOSE-EINSTEIN CONDENSATION OF 
POLARITONS 



In the previous section it was shown that X-P scatter- 
ing is efficient enough to create a quasi-thermal popula- 
tion into the LPB. We will quantify now in detail how the 
occupancy of the lowest energy polariton state evolves 
when the polariton density is increased. 

Above a critical pump-power, the continuous Bose- 
Einstein distribution shown in Fig. || saturates, and 
the additional excitons added by the pump fall into the 
bottom of the LPB. Thus, the lower polariton distri- 
bution becomes discontinuous, and an abrupt increase 
is achieved in the number of polaritons in the lowest 
microscopic level, A^o- Above the threshold for BEC, 
the macroscopic occupancy of a single microscopic mode 
opens the possibility for spontaneous symmetry breaking 
and the appearance of coherence in that mode. The in- 
teraction between the Bose-Einstein condensate and the 
low energy polariton modes would lead then to the re- 
construction of the polariton spectrum, as predicted in 
the well known Bogoliubov approach ||3^. These effects 
can modify the relaxation dynamics above threshold, so 
that our results are not rigorous in this range. However 
we can predict the critical density and temperature for 
BEC, and the most favorable conditions for its observa- 
tion. 

Fig. |l^ (left) shows the evolution of A^o, as a func- 
tion of Px , for different lattice temperatures and polariton 
splittings. For large px above the critical pump, the rela- 
tion between A'o and px becomes linear, and occupation 
numbers of the order of 5 x 10"' polaritons can be reached 
in the ground state. In order to interpret the evolution 
of the ground state occupancy, it is interesting to plot 
A^o (Fig ~ 
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FIG. 10. Left: Evolution of the occupancy of the lowest 
microscopic polariton level, as a function of the pump-power. 
For each polariton splitting we consider Tl = 4, 10, 20/C. 
Right: Same as above, but now the occupancy of the lowest 
microscopic level is plotted as a function of the steady-state 
exciton-polariton density. 

The abrupt increase in Aq allows us to estimate a crit- 
ical density for BEC, which results to be around 
O.lrisQt. However, the situation is now not so clear as 
in the case of equilibrium BEC, where one expects that 
A'o follows the relation (see Appendix B): 



Ao 
S 



(25) 



10 (right)) as a function of the steady-state 



exciton-polariton density, n^^ 



That is, for densities above n'i^, all the excitons added 
to the system should fall into the ground state. Fig. |10| 
(right) shows that, for large densities, the linear relation 
is satisfied. However, close to A^o shows a superlinear 
dependence on Ux-p. This fact can be interpreted as fol- 
lows. The equilibrium BEC relation between A^o and 
would be given by the extrapolation of the linear regime 
of the curves in Fig. |l^ (right), towards A^o ~ 0. By this 
method, we would obtain the equilibrium n^p, which is 
smaller than the one we get from our complete nonequi- 
librium calculation. This means that in a real microcav- 
ity, njp is larger than the one predicted by equilibrium 
BEC. Obviously, this is due to the radiative losses in the 
LPB, and also to the fact that the system has to evolve 
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towards the ground state, following the steps described 
in the previous section. Thus, the X-P scattering has 
to overcome the radiative losses, and the exciton reser- 
voir has to be heated, so that relaxation into the ground 
state becomes efficient enough. The nonequilibrium crit- 
ical density has to be larger than the equilibrium one, for 
these conditions to be fulfilled. 

The transition from superlinear to linear dependence 
of iVo as a function of Px, and n^p, is reflected in the 
evolution of the exciton reservoir temperature, T^. Fig. 

shows that the heating of the exciton gas slows down 
for the pump-power for which the linear growth of A^o 
starts (see Fig. left panel, ftp = lOmeV). The same 
behavior in the evolution of has been obtained for all 
the Op's considered. 

We can also see in Fig. |l^, that for small polariton 
splittings, smaller pump-powers and polariton densities 
are needed for BEC, because the energy distance between 
the bottom of the LPB and the exciton reservoir is also 
smaller. For large Op, the exciton gas reaches higher 
temperatures in the steady-state regime, and the X-P 
scattering rates are slower (thus, larger polariton den- 
sities are needed in order to get BEC). However, in all 
cases considered in both figures, the critical density is 
well below Usat- 

Up to now we have assumed Lc = SO/xm, of the order 
of typical spot diameters. However, Lc is not a well de- 
termined quantity in experiments. Equilibrium BEC in 
2-D systems predicts that both critical density and tem- 
perature depends only on the order of magnitude of Lc- 
In order to check what happens in our nonequilibrium 
situation, we have calculated the density of polaritons in 
the ground state, Nq/S, for different quantization lengths 
(see Fig. 0). 




0.05 0.10 0.15 0.20 0.25 

FIG. 11. Evolution of the density of polaritons in 
the lowest microscopic state, as the exciton-polariton 
density is increased, for different quantization lengths, 
Lc = 25, 50 - 125^m, and ^Ip = lOmeV , Tl = lOA'. 

The occupancy of the lowest polariton mode has a 



more abrupt increase for smaller Lc- However n^^ is 
of the order of OAusat for all the quantization lengths 
considered. 

V. CONCLUSIONS 

Our main conclusions are: (i) At large densities, relax- 
ation after continuous non-resonant excitation can be de- 
scribed by a model in which lower polaritons are injected 
from a thermalized exciton reservoir by X-P scattering. 
(ii) Our model shows a linear relation between the energy 
of the maximum in the lower polariton distribution and 
the exciton temperature or pump-power. By increasing 
the exciton density, the heating of the exciton reservoir 
allows the excitons to overcome the phase-space restric- 
tions and fall into the LPB, that becomes a polariton 
trap in reciprocal space, (in) Bose-Einstein Condensa- 
tion is possible for densities well below Ugat for CdTe mi- 
crocavities. (iv) The most favorable conditions for BEC 
are small polariton splittings, due to the smaller differ- 
ence between the bare exciton energy and the bottom of 
the LPB. (v) Our results have a weak dependence on the 
quantization length, which is taken as the excitation spot 
diameter. 

Up to now, non-resonant experiments performed in 
GaAs microcavities confirm that the weak coupling 
regime is reached before X-P scattering allows polaritons 
to relax towards the ground state |l^Jl^. Experimental 
evidence for the intermediate stage in the evolution of 
the polariton distribution has been achieved in a GaAs 
microcavity pO| ] , and the experimental results agree very 
well with the results presented in section 111. In CdTe, 
experimental evidence for stimulation towards the lower 
polariton states has been achieved under pulsed non- 
resonant excitation. The angle-resolved measurements 
after pulsed excitation in CdTe presented in quali- 
tatively agree with our results for the evolution of the 
polariton distribution as pump-power is increased. How- 
ever, the physical situation in pulsed experiments is quite 
different than the one considered here, because the sys- 
tem is not allowed to reach the steady-state in the re- 
laxation process. Continuous wave excitation seems to 
be a cleaner experimental situation for the study of the 
evolution of the polaritons towards a quasi-equilibrium 
distribution. Our model allows to predict not only the 
possibility of BEC, but also quantitatively describes the 
evolution of the system towards the ground state. Angle- 
resolved measurements under continuous non-resonant 
excitation would allow to show the macroscopic occu- 
pancy of the k — mode, and also to monitor the evolu- 
tion towards BEC as pump-power is increased. 
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APPENDIX A: CALCULATION OF THE X-P 
SCATTERING RATES 

The calculation of the rates for the scattering of ther- 
malized excitons into the LPB was performed in [^ . 
Here we reproduce, for clarity, the main steps, and cal- 
culate also the rate for ionization out of the LPB. 

We consider first the process that injects polaritons. It 
is given by: 



4 



^^31 — K I fx I : -J 1 t^, tj, ; 



(Al) 



iVf , iVJ are the occupation numbers for the initial exci- 
tons in the reservoir, and A^^, correspond to the fi- 
nal exciton and lower polariton, respectively. The in 
the thermalizcd exciton reservoir follow the relation (|^). 
Due to their large density of sates we can neglect the 
stimulation terms in the exciton reservoir (JV^ << 1). 
Energy conservation is implicit in Eq. (Al), so that 
tf ^ ei + ej — . The function R describes the phase- 
space restrictions for the problem of exciton-exciton scat- 
tering in 2-D, under the approximation of isotropic ex- 
citon distributions. It can be expressed in terms of the 
corresponding in-plane momenta |18| : 



^ ■ ' [(fcf+fc'P)2-x] [x-(fcf-fc'P)2] X 



-dx I 



f) 



-x\ \x- 



-{k^^-kff 



(A2) 



where integration is limited to the interval where the ar- 
gument of the square root is positive. A great simplifica- 
tion is obtained when we take the limit /c'^ « k^: 



1 



(A3) 



If we insert this expression into (Al), we obtain 

271 



Wl 



hiksTx) 



1 



which after integration in energy gives the expression re- 
ported in Eq. (ph. 



The calculation of the rates for scattering out of the 
LPB follows the same way: 



m'pl I R{e),e'^, ef , ef)NfNl^de^de^. (A5) 



The factor R is again given by ( |A3| ) , because phase-space 
restrictions are the same as in the previous case, 



Px 



h-KikBTxY 



ej can be expressed in terms of ef , e|, and the expression 
in Eq. ( p^ is obtained. 



APPENDIX B: BOSE-EINSTEIN 
CONDENSATION IN TWO DIMENSIONS 

Let us consider an ensemble of free bosons that fol- 
low a Bose-Einstein distribution with chemical potential 
/i < 0, and finite temperature T. The density of particles 
can be evaluated as: 



1 f°° 
n{p,T) = ^ 



de 



g(£-M)/feBT _ l 



(Bl) 



In 3-D, p(e) cx -^/e, and the continuous Bose-Einstein 
distribution with /i = gives a finite critical density 
nc{T). For densities n > nc{T), all the extra particles 
fall into the ground state. 

In 2-D, p = m/fi (constant), and the integration in 
Eq. (Bl) diverges as /i 0. For an infinite system, BEC 
is thus not possible at finite temperature. For a finite 
system, we can introduce a quantization length, ic, such 
that the integration in Eq. (Bl) has a low energy cut off. 



i{p.,T) = 



27rn^ 



dt 



(B2) 



For the case ec/ksT << 1, Eq. ( |B2D gives a very simple 
result for the maximum density that the 2-D B-E distri- 
bution can accommodate: 



rnksT I e, 
UciT) = — tTT'-°9 



27rr 



knT 



(B3) 



Whenever condition n > nc{T) is satisfied, the extra par- 
ticles added to the system fall into the ground state: 



No{T) 



n,{T). 



(B4) 
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For free particles, ec cx: ^, and Eq. (B2) allows to pre- 



dict a weak logarithmic dependence of the critical density 
on the quantization length. This result qualitatively ex- 
plains the weak dependence of the critical density on Lc, 
obtained in the case of nonequilibrium BEC discussed 
above. 
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